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Method of derivation of phenomenological thermodynamic potential of solid solutions is proposed 
in which the interaction of the order parameters of constituents is introduced through the account of 
elastic strain due to misfit of the lattice parameters of the end-members. The validity of the method 
is demonstrated for Ba^Sri-^TiOs system being a typical example of ferroelectric solid solution. Its 
phase diagram is determined using experimental data for the coefficients in the phenomenological 
potentials of SrTiOs and BaTi03. In the phase diagram of the Ba^Sri-^TiOa system for small Ba 
concentration, there are a tricritical point and two multiphase points one of which is associated with 
up to 6 possible phases. 
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I. INTRODUCTION 

The solid solution of BaTi03 with SrTi03 is one of the 
best documented systems of ferroelectric solid solutions 
(see Refs. [1-3] and references therein). 

BaTiC>3 (BT in short) is well known as the first oxygen- 
octahedra ferroelectric of the perovskite structure with 
paraelectric cubic phase Pm3m and three ferroelectric 
phases: tetragonal P4mm, orthorhombic Amm2, and 
rhombohedral R3m. All three transitions are described 
by one three-component order parameter (polarization), 
and with the potential expansion up to the sixth order 
[4] . Best agreement with the experiment is achieved when 
not only one constant in the term quadratic in the order 
parameter is assumed to be temperature-dependent [5]. 

SrTiOs (ST in short) represents a textbook example 
of a "quantum paraelectric" [6] with a polar soft mode, 
which has never condensed even at very low tempera- 
tures. At a temperature T a of about 105 K pure ST 
undergoes an antiferrodistortive (improper ferroelastic) 
phase transition Pm3m — » I4/mcm with doubling of the 
primitive cell volume. Microscopic mechanism of the 
transition is attributed to the instability of a soft lattice 
mode associated with out-of-phase rotations of TiOg oc- 
tahedra. The transition is phenomenologically described 
by three-component order parameter transforming ac- 
cording to i?25 representation from the point R ( 1/2, 
1/2, 1/2) of the Brillouin zone (BZ) of a simple cubic 
lattice with the Landau potential of the fourth order [7] . 
The constants of the potential were found in Refs. [8-10] . 

Along with softening the soft mode from the BZ 
boundary on cooling [11, 12], the static dielectric permit- 
tivity significantly increases [6, 13] due to softening of IR- 
active F\ u mode from the BZ center [14, 15]. Down to 1.5 
K, the ST crystals remain in the tetragonal phase [6, 16]. 
To describe the temperature behavior of the static dielec- 
tric permittivity, the Barrett quantum equation [17] is 
applied [6]. Thus, the phenomenological potential, which 
simultaneously describes the above-mentioned properties 



should include two three-component order parameters 
-R25 © F\ u . The constants of the fourth order poten- 
tial with such order parameters were detemined in Refs. 
[18-20]. 

According to a Curie- Weiss law in ST, wich follows 
from dielectric and optical measurements [6, 14, 15, 21, 
22] , a ferroelectric phase transition should occur at a tem- 
perature near 35 K associated with softening of the IR- 
active mode. No dielectric anomaly at this temperature 
was found but many authors observed elastic anomalies 
[23-30]. However, in careful dilatometric studies, no low- 
temperature anomalies were observed [31]. 

The solid solution Ba^Sri^TiOa (BST in short) of 
these two remarkable compounds attracts attention of 
researches from at least two aspects: 1) basic interest 
concerning such problems as impurity induced ferroelec- 
tric phase transition in the ST incipient ferroelectric, as 
interplay between the antiferrodistortive and ferroelectric 
phase transitions, as development of the theory to pre- 
dict the (x — T) phase diagram; 2) applied aspect as the 
possibility to control properties and the phase transition 
temperature of these promising for various application 
materials. 

At first, the Ba^Sri-^TiOs solid solutions have been 
studied at 1 > 0.1 and at temperatures more than 
about 100 K. Hegenbarth [32] and Bednorz [33] stud- 
ied imputity-induced ferroelectric phase transitions in the 
BST solid solutions at small x and down to liquid helium 
temperatures. Bednorz [33] and Miura et al. [34] inves- 
tigated the effect of Ba substitution on the antiferrodis- 
tortive phase transition. Later on, this system at low x 
was studied in more details [35-38]. 

The summary of all the experimental data concerning 
the (x — T) phase diagram of BST is the following. 

The temperatures of all three ferroelectric phase tran- 
sition decrease with x decreasing with different rates and 
probably converge at x around 0.1. 

At x = 0.2 + 0.25 [35] on = 0.4 [39], there is a tri- 
critical point for the cubic - tetragonal phase transition 
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when the first order phase transition with x decreasing 
transforms to the second order one. 

The antiferrodistortive transition temperature T a de- 
creases with x increasing, and this transition disappears 
when T a should become lower than ferroelectric transi- 
tion temperature [38]. 

At small x(x = 0.005 [33, 36] and x = 0.01^-0.03 [32, 
35, 38] ), the BST system becomes ferroelectric, and the 
transition temperature T c increases with x increasing. 

In spite of this rich experimental data on the BST 
(x — T) phase diagram, no theoretical treatment of the 
phase diagram is available. This problem is, as men- 
tioned above, of great interest both basic and applied, 
and in this work we develop the phenomenological the- 
ory and derive the concentration phase diagram for BST 
solid solutions. The region of small concentration of Ba 
(small x) is especially interesting since just in the re- 
gion of small concentration: i) the impurity-induced fer- 
roelectricity occurs; ii) there is a strong interplay between 
ferroelastic (antiferrodistortive) and ferroelectric phase 
transitions; iii) experimental data on the symmetry of 
low-temperature phases are not available. 

The paper is organized as follows. After Introduction, 
where background of the problem is depicted, Sec. II fol- 
lows with derivation of the phenomenological potential 
of the BST solid solution. In Sec. Ill , the (x — T) phase 
diagram of the BST system is analyzed, and specific fea- 
tures of the phase diagram are discussed in Sec. IV. The 
last Section presents a summary of the work. 



II. PHENOMENOLOGICAL POTENTIAL OF 
THE Ba x Sn_ x Ti0 3 SOLID SOLUTION 

A. Formalism 

Thermodynamic description of solid solutions we shall 
make by the expansion of the Helmholtz free energy in 
powers of order parameter of a high-symmetry phase [40] . 
At all temperatures, we treat the system as a homoge- 
neous solid solution, which does not exhibit the ordering 
of its components. The constants of the free energy ex- 
pansion are temperature- and concentration-dependent. 
We suppose that we know the potentials of the solid so- 
luton end- members with concentration x = and x = 1. 
Moreover, we suppose that at least at high temperatures 
the end-members exist in one and the same polymorphic 
modification. We shall derive the thermodynamic poten- 
tial of the solid solution of arbitrary concentration using 
the following approaches. 

Atom substitution with no phase transition leads to a 
homogeneous change of crystal size. We shall treat the 
crystal of solid solution as a mixture of two interacting 
crystals with concentration x = and x = 1 with molar 
contribution (1 — x) and x, respectively. The interaction 
of such 'crystals' will be described not by usual introduc- 
tion of additive terms in the thermodynamic potential 
but through the account of elastic strains remaining in 



the framework of the phenomenological theory. These 
strains for the solid solution end-members u and u\ will 
be determined using two conditions. The first one is the 
fit of the lattice parameter of the solid solution a x and 
that of "deformed" end- members (x = and x = 1). 

Let us express the (inequilibrium) lattice parameter of 
the solid solution a via its equilibrium value a x and the 
volume strain u: 



a = a x (l + u) 

Then, the lattice parameters of the solid solution end 
members are 



ao = a x (l + wo), 
ai = a x (l + ui). 



(1) 



The elastic energy of the crystals x = and x — 1 are 
F and F\ respectively. The strain in the F potential 
may be expressed via the strain u of the solid solution 



a-a a x (l + u) - a x (l + u ) 



u 



a 



a x (l + Mo) 



1 



u 



U — Uq 



Similar relation may be written for the F\ potential. 
As a result, the part of the clastic energy related to the 
crystals x = and x = 1 will be (1 — x)Fq(u — uo) 
and xFi(u — ui) respectively. Since the interaction is 
introduced through deformation, the elastic energy of the 
solid solution can be written as 



(1 - x)F (u - u ) + xFi(u - ui) 

Consequently, the second condition is the condition of 
equilibrium obtained by minimization of this energy with 
respect to the strain u. As a — a x in the equilibrium 
(u = 0) we get the equation 



+ x 



u=0 



du 







(2) 



u=0 



It means that the resulting stress should be zero when 
lattice parameter of the solid solution acquires its equi- 
librium value a x . The equations (1), (2) allow to find the 
values of deformations tto and u\, as well as the equilib- 
rium value of the solid solution lattice parameter a x . 

Writing the Fo and F\ potentials in the solid solution 
lattice by means of deformanion, one can introduce, along 
with the common strain u, and common order parameters 
rj. As a result, the potential of the solid solution written 
through the potentials of its end-members will be given 
as 



F= (1 - x)F (r),u - u Q ) + zF^u - m) (3) 

To proceed further, one should minimize Eq.(3) with 
the account of Eqs.(l), (2). 
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B. Potential of the BST system and x = 1 mean ST and BT, respectively). The 

Helmholtz free energy per unit volume F± (i= ST or BT) 

We apply the approach developed in the previous sub- 1S g lven as 

[18] 

section to the Ba x Sri_ x Ti03 solid solution (now x = 



Fi = ftj + + <£) + {vi + + <4) + P12, (vlvi + + flvi) + (p? + vl + pI) 

+an,i {Pi+P2+Pi) + <*i2,i {pIpI+pIpI+pIpI) + "in. (P1+P2+P3) + at 12 M(P2+Pl) 
+Pi{Pi +Pl) + pUpi +Pl)\ + ^W,pIpIpI - (PiP? + V2P2 + VsPa) ~ 1 12. M (pI + pI) 

+ fUp 2 l+Pi) + VUPI +P2)] - *44,» ( ( P2f3P2P3 + flf3PlP3 + fWlP\P 2 ) + F Uii . 



Here, ip is the rotational order parameter (oxygen ion 
shift in the R 25 mode related to out-of phase rotation 
of Ti0 6 octahedral); p , the polarization (related to po- 
lar shifts of ions in the F\ u mode); F u s , the deforma- 
tion potential presented below. The superscript u in Eq. 



(4) indicates that the coefficients are taken at constant 
strain. We shall further suppose that the elastic strains 
are not the order parameters, i.e., there are no proper 
ferroelastic transitions. Then, the deformation potential 
should be quadratic in strain 



F u ,i = -{cu,i + 2ci2,»)ajT(ui + u 2 + u 3 ) + ci 2 ^{u\u 2 + uiu 3 + u 2 u 3 ) + \cn ti {u\ + u\ + u§) 
+ \c i ^ i {u\ + u\ + u%) - 6n,i(ui^f + u 2 <p\ + u 3 ipl) - bi4,,i{unp 2 <p 3 + u 5 ipnp 3 + u 6 ipnp 2 ) 
-6i2,i[tii(v2 + ¥>§) + Mfl + fl) + u 3(fl + - ffn,i(wiPi + u 2 p\ + u 3 pi) 

~944,i( u 4P2P3 + U5P1P3 + U 6 pip 2 ) - gi2,i[ui(p2 + P3) + U 2 (p\ + p§) + U 3 (p\ + 



where Ckj.i are the elastic moduli of the i-th com- 
ponents; a,, their thermal expansion coefficients (i = 
ST, BT); T, the absolute temperature, u is the strain: 
= 8(4x62 k= i 2 , 3; u 4 = + ^£21, . . . , in the 

Voigt notation . Only linear in strain and quadratic in 
the ip and p interaction terms are included in Eq. (5) 
[18]. 

Now, we write the phenomenological potential of the 
solid solution. From Eq. (3), it follows 



Fbst = (1 - x)F S T(p,P,u k - A S T,u m ) 
+ xF BT (ip,p,u k - A B T,u m ) 



(6) 



where k = 1, 2, 3, m = 4, 5, 6, and Ast, Abt are 
determined by Eqs. (1), (2) where u — A ST , Ui — A BT 
with the deformation potentials Eq.(5) at ip = 0, p = 0. 
In the high-symmetry phase Eqs. (1), (2) yield 



a x = 
Ast 



(l — x)TasT+xa.BT 

(l-x)T+x-l(l-x)T+x(l+f)]as T T 

_ -xS-[(l-x)T+x(l+-f)]a S TT 



(l-x)r+x(l+5) 

L+g)[(l-x)rH 
(1-x)t+x(1+5) 



(7) 



a _ (l-i)i-(l+i)[(l-x)T+x(l+7)]tt ST T 
L±BT — " — 



where 6 = 2lSX - 



cii,st+2ci2,st 

Cll,BT + 2ci2,BT ' 
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and a st, &bt are the thermal expansion coefficients of 
BT and ST respectively. 

At room temperature, 5 = 0.026 [35] and 5 = 0.024 
[36]. Below we use <5 = 0.026. The thermal expansion co- 
efficient determines the S temperature dependence. How- 
ever, in a temperature range of - 500 K the change of S 
is not higher than several percents, so we will neglect this 
dependence and consider the 5 coefficient as independent 
of temperature. 

Note that Eqs. (7) give concentration dependence of 
a which differ from the Vegard rule (1 — x)asT + xa B T 
[35, 36]. However, the ratio of elastic moduli in our case 
is close to 1 (r = 1.13), so at fixed temperature deviations 
from the simple linear dependence are rather small. Here 
we may note that it is quite evident physically that lattice 
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parameter of solid solution should depend on the relative 
elasticity of its constituents, so appearance of r in the 
Eq. (7) seems to be quite reasonable. 

The potential, Eq.(6) is the Helmholtz free energy for 
BST with the end-member ST and BT potentials with 
the order parameters which are given now for the BST 
lattice. The Gibbs potentials = Fi — au (i = ST, 



BT, a is the stress) may be written by substitution of 
strain u for the quantities determined from the equation 
a = From the form of Eqs. (4) and (5), it follows 

that only constants of the second and fourth order will 
be renormalized. Thus, the quadratic part of the Gibbs 
potential of BST will be 



{< 



+ 1 



1T + x% iBT + (bn.BT + 2b 12 , BT - b n ,ST ~ 2b l2 , ST ) ^A^zg±j] } ^ 
Is, + *< BT + (9n,BT + 2g 12 . BT - g n ,ST 2g 12 , ST ) EilAzteg±£l } p 2 



(8) 



Below we will use the Gibbs potentials 



III. PHASE DIAGRAM OF THE Ba x Sn_ 
SYSTEM 



<Ti0 3 



To study the solid solutions in the framework of the 
formalism developed above, it is necessary to determine 
the potential coefficients at the boundary points x = 
and x = 1 for the &st , and $ bt potentials, respectively. 

As a base for $st, we take the potential up to the 
fourth order presented in Ref. [18] and modified in Ref. 
[20]. Taking into account the quantum properties of $st, 
which arc important at low temperatures, we assume that 
constant (3° ST in the if 2 term from Eq.(8) follows the 
equation [41] 



B 



coth 



60.75 



coth 



60.75 



(J/m 5 ) (9) 



with T a = 106 K . 

With experimental data from [18] 3^ ST = 1.69 x 
10 50 J/m 7 and f = 0.69 x 10" n m at T = 2 K, one ob- 
tains B = 1.7 x 10 28 J/m 5 . Note that this value of if is 
associated with the angle of TiO@ rotation which is equal 
to 2°. 

For the $ bt potential of BaTiOs we take the constants 
from Refs. [5, 42]. The potential with these constants 
satisfactorily describes the phase transitions in BaTi03. 
The BT crystal is stable with respect to appearance of 
the ip order parameter at the R (1/2, 1/2, 1/2) point of 
the BZ boundary. Therefore, we add to the potential [42] 
the part quadratic in ip with the constant proportional 
to the square of the appropriate "soft mode" frequencies 
at the R point. The order parameter if at this point 
is associated with oxygen ion shifts (oxygen octahedral 
rotation) only [43] , and Ba ions do not take part in these 
shifts. Thus, the ratio of this frequency squared to the 
constant in the f 2 term should be the same as that in 
the ST crystal at the same temperature. Taking for the 
appropriate frequencies at T — 300 K the values [12, 44] 6 



meV and 16 meV for ST and BT respectively, we obtain 
B a = 3.7 x 10 29 J/m 5 . 

We neglect the other (/3-dependent terms in $bt as 
compared to the (3° BT f 2 term since there is no transition 
with the if order parameter in BT. 

The coefficient in the p 2 term we take from Refs. 
[20, 45], where this value was found from the temper- 
ature dependence of ferroelectric soft mode with account 
of specific behavior of SrTi03 at low temperatures: 



a a a =A 

1,ST 



, . 54 \ , /54 
coth — — coth — 
T V 30 



Jm/C 2 (10) 



According to Ref. [20], A = 4.05 x 10 7 Jm/C 2 . This 
value gives the Pm3m <-> P4mm transition line in the 
T — x phase diagram which does not agree with exper- 
iment (almost straight line in the experiment). If the 
thermal hysteresis is neglected, this transition line should 
be determined by the line of zero coefficient in the p 2 
term in Eq.(8). We assume that at high temperatures 

In such a case coth(a/T) w T/a, 



Qa. a da 

1,ST 1,BT 



dT 



dT 



and the xT term in the first two terms at p 2 in Eq.(8) 
may be neglected, and the transition line becomes close 
to the straight line. Then, one obtains A = 1.8 x 10 7 
Jm/C 2 . However the best agreement with the experi- 
ment is reached with A = 1.5 x 10 7 Jm/C 2 . Just this 
last value we will use below. 

Constants from Eqs. (9), (10), along with the con- 
stants B a , t a and t a determine the behavior of 

11, ST ' 11, ST 12, ST 

ferroelectric modes at temperatures below the transition 
point T a . Since the a° ST constant was changed, it is 
necessary to redetermine the t° , t° 2 ST constants. 

To do this we used the following procedure, similar 
to that of Ref. [45]. Neglecting a small hard modes 
contributon to dielectric constant we have [3] 

1 = Q^ST = Ruj2 

e dp 2 

K being some constant scale factor, which we deter- 
mine from this relation at the temperatures above the 
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FIG. 1: Temperature dependence of the ferroelectric modes 
in SrTiOa. Experimental points are from Ref. [45], solid lines 
are fitting with t° 1ST = -0.7 x 10 29 , t° 2 ST = -0.32 x 10 29 in 
J/C 2 m units and with scale value K — 2a" x jio 2 = 1.65 x 
10 4 Jm/C 2 cm- 2 at T=200 K. 



ST transition temperature T Q . Below T a analogous re- 
lations (with the same K) for the inverse longitudinal 
and transversal dielectric constants define correspond- 
ingly the temperature dependence of the frequencies of 
the soft Aiu and E„ modes resulting from the splitting 
of the original Fi u mode in ST tetragonal phase The for- 
mulae for these split frequencies are given in Ref. [18, 45] 
and they depend on coupling coefficients t^ ST , t^ ST . 

Using experimental data [45] , we recalculated them to 
obtain t£ ST = - 0.7 x 10 29 and tf a ST = -0.32 x 10 29 
in J/C 2 m. To examine the validity of these values, the 
experimental data [45] on the temperature dependence 
of the ferroelectric modes in SrTiC>3 and the dependence 
predicted using our values of the constants are presented 
in Fig.l for the scale factor K — 1.65 x 10 4 Jm/C 2 cm~ 2 . 



We discuss further our values of the constants ST , 
and ot° 2ST - In the experimental phase diagram [35, 38] 
one observes the convergence of the lines of ferroelectrics 
phase transitions and probably disappearance of the 
P4mm and Amm2 phases. This can be described by the 
sign change of the anisotropic coefficient of the fourth- 
order invariant in &bst which is equal to {pt° — 2a° ). 
Since a a - 2a a > 0, the value of (a" - 2a a ) 

12, BT _ 11, BT ' ^ 13 11' 

in $_bst potential should change the sign as a function 
of concentration x (this follows from the global stabil- 
ity of the BT rhombohedral phase ), which is possible if 
a" — 2a a < in <j>?T- Moreover, one observes dis- 

12, ST 11, ST J± ' 

appearance of thermal hysteresis for the Pm3m «-> P4mm 
phase transition at x = 0.2 0.4 [35, 39, 47], i.e., there 
is a tricritical point in (x — T) phase diagram. 

We assume that three phase transition lines converge 
at x = 0.13 and the tricritical point is observed at x = 0.2 
(x = 0.2 -7- 0.25 in [35]). As a result one obtains a° = 
0.33 and a a 12 ST = 0.2 in 10 9 Jm 5 /C 4 . (If the tricritical 
point occurs at x = 0.4 [39, 47], the constants increase 



TABLE I: Constants of the potentials for the ST and BT 
crystals. Constants with superscript u are for the Helmholtz 
free energy Fi — Fi(<p,p,u), constants with superscript a, 
for the Gibbs potential $i = &i((p,p,cr), (i = ST, BT, u is 
the strain, a, the stress). Constants in the quadratic terms: 
ff ST = 1.7 x 10 28 [cot/i,(60.75/T)-cot/i(60.75/106)], ft BT = 



3.7 x 10 29 


in J/m 5 units, 


a a 1ST = 1.5 x 10 


1 L.UL1 1 1 OI j ± J 


cotA(54/30)], < BT = 3.34 x 10 5 (T - 381) in 


Jm/C units. 


Constant 


SrTiOa 


BaTi0 3 


Units 


"11 


1.94 a 





xlO 50 J/m 7 


A, 

"12 


3.97 a 







1 1 


1.69° 







"12 


3.88 a 





x 10 9 Jm s /C 4 


a" 
1 1 


0.63 (2.0 a ) 


-1.48+0.00469T 6 


„ u 
u 12 


0.02 (1.19™) 


-0.0877 6 




a a 
"n 


33 Cl 7 a ) 


-2 04+0 00469T 6 




a CT 

12 


0.2 (1.37 a ) 


0.323 6 




„ u „ u 
"in ' "m 





24.45-0.0552T 6 


x tu JVSX j 


„ U „ CT 

CK, ,„ . (X 

112 ' 112 





4.47 b 




„ u „ u 
rv fy 

123 ' 123 





4.91 6 


x iu J 1 vn 


1 11 


-2.0 (-3.0 a ) 





''12 


0.67 (0.23 a ) 







l ll 


-0.7 (-1.74°) 







b 12 


-0.32 (-0.75°) 







Cll 


3.36 a 


2.06 c 


xlO 11 J/m 3 


C12 


1.07 a 


1.4 C 




C44 


1.27 a 


1.26 c 




6n 


1.25 a 





xlO 30 J/m 5 


bi2 


-2.5 a 







&44 


-2.3 a 







911 


1.25 a 


l.O 6 


x 10 10 Jm/C 2 


312 


-0.108° 


-0.017 6 




044 


0.243 a 


0.365 6 





Ref. [20], 6 Ref. [42], c Ref. [46]. 



up to ofj ST = 0.69 and of a ST = 0.92 in 10 9 Jm 5 /C 4 ). 

All the constants in the <I>st and &bt potentials ac- 
cepted in this work are presented in Table I. 

For the constants presented in Table I, the coordinates 
of the convergence point are given as x = 0.13, T = 93 K, 
and these of the tricritical point are x = 0.2, T = 118 K. 
In the calculated phase diagram there appears one more 
N-phase point with the coordinates x = 0.027, T = 48 
K. The theoretical phase diagram together with experi- 
mental points [35-39] is shown in Fig. 2. The calculations 
both analytical and numerical were performed with the 
Maple-6 software package. 

Symmetry of the phases and their positions near the 
second N-phase point depend on the constant fj 4 which 
seems to be not quite well determined in the experiment 
contrary to t\ x and t\ 2 - Depending on the sign and value 
of £44, the following three phases can appear near the N- 
phase point with temperature decreasing: ((p,ip,tp,p, p, p) 
- R3c, (Q,0,tp,p,p,0) - Ima2, and (ipi,ipi,ip2,Pi,Pi,Pu) ~ 
Cc. 

With negative t% 4 , only the phases Ima2 and Cc exist 
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FIG. 2: Phase diagram of BaxSri-^TiOa Solid and dashed lines present calculated diagram for the first- and second-order 
phase transitions, respectively. Experimental points: circles [35, 38], crosses [37], triangles [36], squares [39]. The calculated 
coordinates of the tricritical and multiphase points are Tfc = 118 K, x — 0.2; Tjvi=93 K, x = 0.13; Tjv2=48 K, x — 0.027. 
Detailed pictures of the phase diagram in the vicinity of the multiphase point Tn2 are presented in Fig. 3. 



(Fig. 3a), and transitions to the Cc phase are only of first 
order. At small positive £44, the phase R3c appears in a 
narrow region between the phases Cc and R3m (Fig. 3b). 
The phase borders between R3m <-> R3c and R3c «-> Cc 
phases are second-order. As £J 4 increases, the region of 
existence of the R3c phase becomes broader because of 
decreasing of the existence region of the low-symmetry 
phase Cc unless the latter disappears completely. In this 
case, the transition line between the R3c and Ima2 phases 
is of first order (Fig. 3c). As £44 increases further, the 
Ima2 phase becomes "swallowed up" by the R3c phase, 
and first-order phase transition is possible between the 
I4/mcm and R3c phases (Fig. 3d). At £J 4 > 3.1 x 10 29 
J/C 2 m, the (x = 0) potential describes the R3c 

phase at low temperatures, which at first is metastable 
and then at £44 > 4.0 x 10 29 J/C 2 m becomes stable. 



IV. DISCUSSION 

Now, we discuss the phase diagrams in Figs. 2, 3. 
One can see that in the region of small Ba concentra- 
tion x, there is rather a large scatter of the experimental 
points. The scatter seems to be due to different meth- 
ods of sample synthesis. E.g., changing temperature of 
the synthesis, one obtains different dielectric properties 
in the samples with the same concentration x [48] which 
should lead to different values of the constant of the po- 
tentials. In Fig. 2, the experimental point at x = 0.1 
from Ref. [37] cannot be used, otherwise a\\ t sT should 
be negative. In this case more complete model is needed 
with up to the sixth order in <E>st- 

In the phase diagrams in Figs. 2, 3, the existence of N- 
phase points T^ri and Tjv2 do not contradict the Gibbs 
phase rule. This follows from the developing of the phe- 
nomenological theory itself: along the lines and in the 
points of second-order phase transitions, a system is in a 
single-phase state, namely in the phase with the highest 
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FIG. 3: Phase diagram of Ba^Sri-^TiOa at low x in the vicin- 
ity of the Tjv2 = 48 K multiphase point. As in Fig. 2, solid 
and dashed lines present calculated diagrams for the first- 
and second-order phase transitions, respectively.Experimental 
points are from Refs.[35, 38] (circles) and from Ref. [36] (tri- 
angles). The values of tj 4 in figures a)-e) are: a) < 0, b) 
< < 2.3, c) 2.3 < C < 2.8, d) 2.8 < £ < 4, e) £ > 4 
(in x 10 29 J/C 2 m units). 



symmetry. One should apply the Gibbs rule only to the 
solid lines in Figs. 2, 3 (first-order transitions); in the 
points of these lines, the system is multi-phase, in our 
case, two-phase. 

Various versions of the phase diagram at small x are 
presented in Fig. 3. The steepness of decreasing the 
1 (Pm3m <-> I4/mcm) phase transition line as a 
function of x is determined by the frequency of the ip - 
normal mode at the R (1/2, 1/2, 1/2) point from the BZ 
boundary which in the BT crystal is high and is assumed 
to be temperature-independent. 

For the line /3\ = (the «-> 1 phase transition line), 
one obtains x — 0.034 at T = 0. If the position of the 
0^2 (Pm3m «-> P4mm) phase transition line (see Fig. 
2) changes depending on the sample synthesis, the N2- 
phase point would change its coordinates along the line 
Pi = 0. If TjV2 decreases, the region of stability of the 
phase 1 (I4/mcm) will increase up to x = 0.034. If T/V2 
increases, there is a possibility of crossing the borderline 



of the phases 5, 6, 7 and the line x — 0. Thus, several 
phase transitions in the pure ST crystals become pos- 
sible. Impurities change the deformation potential [49] 
and that results in renormalization of the constants of 
the complete potential and in the shift of points in the 
phase diagram. 

Let us pay attention to the phase borderlines in Figs. 
2, 3 which are close to vertical. Even with a very precise 
determination of the Ba concentration, there is always 
concentration ingomogeneity. That is why the vertical 
borderlines will be smeared. We would like also to note 
that another disturbance of the phase diagram may be 
due to possible spinodal decomposition of the solid so- 
lution and formation of BT nanoregions in the process 
of sample preparation [50]. Note that maybe just such 
a decomposition of the solid solution was the reason of 
glasslike behavior in the ceramic samples at low x ob- 
served in Ref. [35] , though this result was not supported 
in the other works [33, 36, 38]. 

Now, we discuss how important is the elastic contribu- 
tion for the phase diagram. It appears that at, for ex- 
ample, x = 0.51, the difference between the 0-2 phase 
transition temperatures calculated with and without the 
elastic contribution is 24 K. At x = 0.8, the same value 
for the 3-4 phase transition temperatures is 17 K. ( The 
elastic contribution decreases the phase transition tem- 
peratures.) Thus, the difference is no more than 10%, 
i.e., not very high. However, even small variation of the 
potential constants can significantly change the phase di- 
agram at small x, what is clearly seen in Fig. 3. 



SUMMARY 



The method of construction of the potential in the 
framework of the phcnomenological theory has been de- 
veloped for description of the concentration phase dia- 
gram of non-ordering solid solutions. ( The ordered solid 
solution is a special case which should be treated sepa- 
rately.) 

The concentration [x — T) phase diagram of 
Ba^Sri-^TiOs solid solution has been derived, and the 
diagram adequately agrees with the available experimen- 
tal data. 

There are three specific points in the phase diagram: 
two N-phase points and one tricritical point. The coordi- 
nates of these points critically depend on the numerical 
values of the potential constants, which may explain large 
scattering of experimental data. 

At low Ba concentration, there are several phase states, 
which are very close with respect to the concentrations. 
Concentration fluctuations in the samples can result in 
broadening of these " quasi- vertical" borderlines. 
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